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Introduction and sculpting the problem. 

In this work attempt to generalise Bogoliubov's chain, 0, on relativistic case, when the particles 
number is not fixed, proceeding from the QFT dynamic equations is perfomed. Our consideration will 
be based on the Wigner's function formalism. 

1. Notations and preliminaries. 

In d-dimensional quantum mechanics, if the system is in a pure state with the wave function the 
Wigner's function looks as 

F(k, x) = -. 1 w / dve* +kv V*(x + -v)V(x - -v). 
K ' ; {2nh) d J v 2 J K 2 J 

In relativistic quantum statistical physics this formula transformes as ,0: 

T{K X) = (2^p / dVe ^ V < + 1 V)IP{X ~\ V)> ' (L1) 

where (p* and tp are the secondary-quantizated field operators and averaging is performed by Yuttner- 
Singh distribution, see below. According to (1.1) the one-particle Wigner's function is the average of 
filling number, then the norm condition is: 

J J r ( y k,x)dkdx = N (1.2) 

Let's note through z% = (ki,Xi) the coordinate of i-th particle in relativistic quantum-mechanical 
phase space. 

Let's inductionally introduce n-particles Wigner's function in relativistic case: 

Tj^x . . .z n ) = W J exp fi^kiVij W 2n (xx + y , . . . ,x n - y^ dvi . . .dv n . (1.3) 
as the Furier-transform of the normal-ordered and averaged operators product: 

VV2n [Xl + —,... ,X n —I — 

='<■■ (*. + ^)..V (*. + £)*>(*!-£)...*>(*.-£):>, (1-4) 

were 

N' 1 



(2n) dn ' 

Heretheto we imply h = c = k^ = 1 and field operator tp(x) is taken in Haizenberg representation: 

<p{x) = U(tYT(cp (x)U(t)) (1.5) 

,where <po(x) is the solution of the uniform Klein-Gordon's equation and U(t) is the evolution operator: 

t 

U = Texp(i J H int (t)dt) (1.6) 
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Averaging-out in (1.2) is produced as to ordinary rules of statistical mechanics: 

<: O :>= Tr(p: O :). (1.7) 

Here trace is calculated as to any full set of states. The statistical operator p is the relativistic analog 
of the operator realisation of the grand canonical Gibbs distribution, where the fact, that the particles 
can produced and disappear, is allowed for by Haizenberg's representation, see(1.3): 

pit) = ±exp((3(pN(t) - U v P"{t)), (1.8) 

where N(t) is operator of number of particles in Haizenberg's representation: 

N(t) = f d 3 p[ai(p,t)a(p,t) +tf(p,t)b(p,t)\, (1.9) 



where (3 = ^, p is the chemical potential; U v is the hydrodynamical velocity 0, and ( is defined by the 
norm condition: Tr(p(t)) = 1. 

The 4-momentum operator of the system is 



,where k = 1, 2, 3 and 



P»(t) = (H(t),P k (t)) 



H(t) = / d 3 pp°[a\p, t)a(p, t) + b\p, t)b(p, t)} (1.10) 



, where p° = ^/p 1 + m 2 . Here and are operators of particles and antiparticles production in 
Haizenberg's representation correspondingly, see(1.5). 

Amenably (1.2) the norm condition for s- particle Wigner's function takes form: 

J Fsifa, x 1 ...x s , k s )dz 1 ...dz s = N s , (1-H) 

where iV is the general particles number in the system. 
2. The model description. 

In this work we'll consider the quantum statistical relativistic system, consisting from elementary 
particles, interacting by Lagrangian of self-interacting complex scalar field: 

C = (<9^<9^)m%V) + A(V^V) 2 (2-1) 
On the strength of corresponding this model field equations: 

(□ + m 2 V(x) =p(x); (2.2) 

(□ + m 2 )^{x) = p\x) (2.3) 

, where 

p(x) = 2\(<p\xMx))<p(x); (2.4) 

p ] = 2\(ip\x)ip(x))ip j (x). (2.5) 
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we'll deduce BBGKY hierarchy below. 

3. The kinetic equation for n-partical Wigner's function. 

The Lorentz-invariant kinetic equation for one-particle Wigner's function, as to essence, was already 
obtained in the monography[[|, and look as: 

*mW) = / elkV < P( x + lM x ~ \)) > dv + h - c - (3- 1 )- 

By substitution like-current field p (2.5) in the given equation we have: 

k^F{z) = J <■ vH x + gM* - 2 ))¥? (x + 2 + 2^ X + 2 ~ 2" )} :> X 

Xe -ikV e -ikv dk > dvdv , + hc (3.2) 
Here we took into account the integral representation for 5-function: 

So, we have received the first equation of engaging equations: 

k^d^iz) =i\Jf 2 (k, x, k', x + 7-) e~ lk ' v dk'dv + h.c. (3.3) 

The subsequent equations of chain could be found by induction, i.e. it is valid the next 
Proposition 1 

kf J ,d fl \ Xi J' a (zi, . . . z n ) 
= ixj T n+1 {z x , ...z n , z')e- ik ' v dvdk' + h.c, (3.4) 

where z' = Xi — |v, k' . 

4. The renormalization of the chain equations. 

As we well know, in QFT in the force of it's locality singularity arise at operator-meaning generalized 
functions(field operators) production. For their eliminating one resort to renormalization trick in order 
to radiation corrections have not led to infinite expressions for the couple constant and mass, then in 
order to search renormalization of the received chain it is essentially to answer on the following question: 
"How to distinguish renormalized (observed) part of the n-particleal Wigner's function at transition to 
the semi-QFT (relativistic thermodynamics)?" 

For answer on this question advisably to express an arbitrary operator in terms of reduced density 
matrices describing initial state, then the Hilbert space of the system is defined by initial states, which 
are identical with asymptotical ones (t — > — oo for infinity remote past). 

For the first we shall consider, for simplicity, the case without antiparticles, then the asimptotic 
pure state is: 

\p n ) hl = ai(p n )\0), n = 0,l,2,..., (4.1) 

Here we have used the compact notification: p n for set Pi,P2, ■ ■ ■ ,P% and a\ n {p n ) for asimptotic 
operators product a\ n (pi) . . . a\ n (p n ). 
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If we cut-off Hamilton operator by volume and tribute the periodic boundary conditions then we 
have the discrete energetic spectrum corresponding to Hamiltonian: 



H(t)= / d 3 pp°[J(p,t)a(p,t)+tf(p,t)b(p,t)). 



Pmin 



The number of particles is finite and is eigenvalue of operator 



Pmax 



N(t)= / d 3 p[a\p,t)a(p,t) + b\p,t)b(p,t)]. 



According to that the volume of our system is very small quantum effects play essential role. Then 
energy of system as due to Haizenberg's ambiguity principle is very large and upper integration limit as 
to momentum p max could be entered as particle momentum, generated after interaction immediately. 
The lower integration limit p m i n could be entered as particle momentum in the time moment of system 
evolution accordihg to Haizenberg's ambiguity principle: E ~ -, where r is the life-living time of the 
emmited secondary particles. We stress here, that in this situation statistical considerations are applied, 
i.e. the number of particles is large, but finite, unlike QFT for which, as we well know, infinitely large 
degrees of the freedom (field quants) is characteristic. In force that the number of particles is finite the 
completeness condition is next: 

N{E) 1 d 3 n 

1= E^/J^b n >*in<n (4-2) 



n=0 



n\ J p 



where N(E) is the function of energy of the system. Then the modification of formula for average for 
arbitrary operator O, ref.[§], with take into account new completeness condition(4.3)is neaded: 

N(E) . N(E) / i\m r n+m 

(O) = Y, ~, / d A x n d A k n O n (x n ,k n ) ]T / d 4 x m d A k m J] ^{x^kj). (4.3) 

n=0 n - J m=0 \ m -) J j=l 

Here the integrals as to momenta are cutted off on upper and lower limits according with above said 
and integrals as to space are cutted off in accordance with the one meson space volume. Multiplier 

N(E) m 

£ t^tyt J d 4 k m d 4 x m n F- m {xj,kj) equals 

m=l ^ '* j=l 

Z = -< N m (E) > +\(< N m (E) >) 2 - i(< N in (E) >) 3 + . . . + fe^(< N m (E) >) N{E) , (4.4) 
where we used, that 

J Ti n {k, x)dxdk =< N in > (4.5) 

. Ni n is the number of scalar particles operator in initial state. Let's examine the Wigner's function 
expansion as to the contributions of collisions. After substitution (4.5) into (4.4) we have: 

n(e) 1 

(0> = E -7 / d A x n d A k n O n (x n ) k n ) J] F in ( Xj , kj) x Z (4.6) 

n=0 n - J j=l 
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Having applicated (4.6) to the one-particle Wigner's function, we have the following expansion: 

n(e) 1 

^(x,k) = E - / dVdV^n(x I1 ,p I1 |k,x)n^in(x + x J ,p J )Z (4.7) 

n=0 n - J j=l 

, where 

r n 1 1 

* n (x n ,p n \k) =N'J d% n J] e^ x \ n {p n - -u n \V(k)\p n + -u n ) m . (4.8) 
,and operator \I>(£;) is: 

*'*»=(2^/^ ,lv G")K4")- < 4 - 9 > 

Let us calculate every matrix element fy n (x n ,p n \k) separately. Matrix element equals zero at n — 0, 
as it have to be for normal-ordered operator. 

We have used the completeness condition (3.7) at n — 1, which leads to the following expression: 

X I d*uexp(rux) ia (p - i<V(0)|p"")) o „ t „„,<p"X0)|j> + iu)^ (4.10) 
Matrix elements — |-u|(/?'''(0)|p" n ) ou t can be calculated by using Jang-Feldman equation. 

Finally we have 

^ 1 (x,p\k) = 5 {A \k-p)5 {A \x). (4.12) 

Proposition 2 

Using the results (4. 10)- (4. 12), the expansion of Wigner's function as to the contributions of collisions 

is: 

N(E) 

^ /dVdV^n(x n , P n |k)x 

n=2 



N(E) 1 

^(x,k) = Z x [^ in (x,k) dVd 4 p n ^n(x n ,p n |k): 



n 

xn^ +) (x + x J)Pj )]. (4.13) 
j'=i 

It is impossible to overlook obvious similarity received expression with virial expansion for the classical 
distribution function in nonrelativistic statistical physics. 

Remark: According to above mentioned assumption about initial state of the system, we see, that 
the first term corresponds to solution of uniform or colisionless kinetic equation: 

k^F m (z) = 0. (4.14) 
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Initial distribution function is the equilibrium relativistic Boze-Einstein distribution function as it dis- 
cussed above, but without spin: 



The rest of the members allow for the contribution of interparticles interactions. According to (4.5) Z 

is the finite value unlike QFT, where it would be diverging quantity and accords to renormalization of 

_i 

wave function: (p R (x) = Z v 2 (p(x), see0 due to Dyson's transformation. Thus the following statement 
is valid: 

Proposition 3: 



Fn\Z\ . . . Z n ) 



iN'ZX I F n+1 {z l ...z n ,z')e- lkv dvdk' + h.c. (4.16) 



' nun 



5. The physical illustration. 

Received chain may be applicated for the computing of kinetic coefficients in non-linear approach, i.e. 
when the effect of correlations between particles is essential. By other words,it is necessary to allow for 
not only interaction of particles with each other, but also the interaction of particles with clusters. For 
example, such situation could be emerged in description of phase transition in relativistic plasma. The 
role of clusters play numerous short-living resonances in the elementary particles physics. Conjecturally 
adequate mathematical tool for non-equilibrium interactions between clusters (resonances) description is 
the Bogoliubov's chains method(BBGKY). Let's for definity consider the following many body processes: 



PP — > pp7Y + -K + 7Y 7T 7T° — > pA ++ 7T + 7T 7T 7T° 
PP — > pnn + 7T + 7T + 7T~ 7T~ 7T~ —> pA + 7T + 7T + 7T~ 7T~ 7T~ (5-1) 

and etc. The volume, within that collision occurs and the number of generating particles could be 
calculated, for example, from multiperipheral model, [||]. In this article it was built the quantitative 
theory of the multiple mesons production with respect of secondary particles interactions at ultra-high 
energies. There it was obtained, what the number of particles logarithmicaly depends on energy E 
generated by collision: N(E) ~ Ins, where E ~ y/s. Respected that volume of system is Lorentz- 
contracted V ~ ^3"^; where M is nucleon mass, maximal tt- meson energy-momentum equals 

E Ins 



Pmax = (Pmax = \l fraax + m 2 ~ — ~ "-^,p max ). (5.2) 

In the case of nucleon beams collision all emitted energy is exuded in the small volume and spent to 
pionic stars formation and belated to the binded states formation or resonances, see. (5.1). With the 
flow of time the system expands and becomes more dilute. I.e. we have a number of, step by step 
alternating each other descriptions: local or QFT-description, then thermodynamical follows, which 
replies to nucleon plasma on the early stages of it's formation corresponding to the high energies, 
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no more than lOGeV or hydrodynamical description M for ultra-high energies, then kinetic stage 
follows, at which correlations accord essential influence on the system dynamics and eventually passes 
in equilibrium state, corresponding to Yuttner-Singh distribution. 
Summary. 

The kinetic stage of the evolution of the system of scalar particles interacting as X(^Lp) 2 , is described 
by the chain of equations (4.16). 
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